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Abstract 

A nonlinear parabolic equation of the fourth order is analyzed. The equation is characterized 
by a mobility coefficient that degenerates at 0. Existence of at least one weak solution is proved 
by using a regularization procedure and deducing suitable a-priori estimates. If a viscosity term 
is added and additional conditions on the nonlinear terms are assumed, then it is proved that any 
weak solution becomes instantaneously strictly positive. This in particular implies uniqueness for 
strictly positive times and further time-regularization properties. The long-time behavior of the 
problem is also investigated and the existence of trajectory attractors and, under more restrictive 
conditions, of strong global attractors is shown. 
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1 Introduction 

This paper is devoted to the analysis of the following class of fourth order parabolic equations: 

Ut ~ div(b(u)Vw) = 0, (1.1) 
w = Sut - Au + f(u) + 7(u) - g, (1.2) 

on £1 x (0,+oo), ft being a bounded smooth subset of R d , d 6 {2,3}, coupled with the initial and 
boundary conditions 

u\ t =o = it , in fi, (1.3) 
d n u = b(u)Vw ■ n = 0, on dQ. (1.4) 

The function b(u) represents a solution-dependent mobility coefficient that possibly degenerates at 
as a power of u (cf. (|2.ip below), while the sum / + 7 stands for the derivative of a configuration 
potential W . In particular, we assume that / is the (dominating) monotone part, with f(u) ~ u~ K for 
some k > 1, and that 7 is a bounded and globally summable perturbation that accounts for possible 
nonconvexity of W. The coefficient S in (|1.2I) is assumed to be nonnegative, with S > describing the 
presence of viscosity effects. Finally, g is a smooth external forcing term. 
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In the two-dimensional case, problem (|l.lD - (ll.4j) can describe the evolution of some classes of 
thin liquid films, with u representing the height of the film. Then, the singular behavior of / near 
accounts for the presence of short-range repulsive forces, while the nonmonotone character of W 
at oo (given by the term 7) is related to the occurrence of long-range repulsive forces. An extensive 
presentation of the underlying physical situation is given in [DJ to which we refer for more details (see 
also [TS] and Remark 12.21 below) . 

In the three-dimensional case, the model is also physically relevant since it is closely related 
with the Cahn-Hilliard equation [8] with nonconstant mobility analyzed in a number of recent papers 
both in the nondegenerate and in the degenerate case (cf., e.g., [3] 2] [T3J HH1 HI] and the references 
therein). Indeed, if ±1 represent the pure states of the order parameter in the Cahn-Hilliard model, 
then we can modify / into the form f(u) ~ (1 — u 2 )~ K and accordingly suppose that b degenerates near 
±1 (instead that near 0) as a power of (1 — u 2 ). Then, the so-modified system (|l.l[) - (ll.2p turns out 
to represent a variant of the Cahn-Hilliard equation with degenerate mobility and singular potentials 
analyzed in the celebrated paper [13 . Correspondingly, all the results proved here for (ll.ll) - (|1.2p also 
apply to the Cahn-Hilliard setting with straighforward modifications in the proofs (to be more precise, 
in the Cahn-Hilliard setting we would even have slightly stronger results since it would no longer be 
necessary to take care of the growth of b at infinity). This also motivates the choice of considering 
also the viscous case S > 0, which is particularly meaningful in the context of Cahn-Hilliard models 
(cf. PU, see also [2D]). 

Initial-boundary value problems related to (|1.1[) - (|1.2I) have been addressed in a number of 
recent contributions. In particular, in [5J further qualitative properties of the solutions are proved 
in the one-dimensional case and the stability properties of the steady states are investigated. The 
papers [15] (devoted to the one-dimensional case) and [14] (considering space dimensions 2 and 3) 
analyze problem (jl.ip - (ll.4[) under assumptions on the nonlinear terms very similar to ours. In [141 115] . 
existence of a solution is proved by means of a nonnegativity-preserving finite-element scheme, which 
is also effective for a numerical investigation for the model. Finally, we mention the recent work 
[23], where the long-time behavior of the problem is studied in the one-dimensional setting. In this 
situation, the authors can prove strict positivity of the solution also in the nonviscous case, which 
allows them to show existence of a smooth global attractor by relying on the standard theory of 
infinite-dimensional dynamical systems. 

Our first purpose in this paper is to prove existence of at least one weak solution to Prob- 
lem (|l.ip - (|1.4p under general assumptions on the data, by using a regularization - a priori estimate - 
passage to the limit procedure. Compared to the proof given in [T3] , our method has the advantage 
to be relatively simple. Moreover, as a byproduct of our procedure we see that the possibly singular 
solutions to (jl.ll) - (jl.4p can be approximated by the smooth and positive solutions of a regular PDE. 
Actually, if we consider, for instance, the nonviscous equation (i.e., the case 5 = 0), then we can con- 
struct the e-approximation taking S £ > (which, as noted above, is physically motivated at least in 
the Cahn-Hilliard setting) and choosing a sufficiently singular function f £ , with 5 e \ and f £ — >• / in 
the limit e \ 0. Then, for e > Theorem 16 . 1 1 applies (cf. also Remark l6.3l) : hence, the approximating 
solutions u £ are smooth and positive. Moreover, the very same argument used to pass to the limit in 
Subsec. I3T31 below shows that any limit point of {u £ } for e \ solves the original nonviscous problem. 

In comparison with |141 115) . we have here the extra assumption that the singularity of l/b(u) 
at u — is not stronger than the singularity of the potential F (the antiderivative of /), cf. (|2.3p 
below. On the one hand, this assumption looks natural and is satisfied for the physically relevant 
examples of b(u) and f(u). Indeed, the most used non-linearity b in the theory of thin films is (cf., 
e.g., © or M) 

b(u) = u 3 + f3 3 - n u n , ne(0,3), (1.5) 
and for the function / we have the following model examples (cf. [5]): 

f{u) = A11- 3 - Bu- 9 or f(u) = Au' 3 - Bu~ A , (1.6) 

where A,B > 0. Thus, in all the physically relevant cases (|2.3p is satisfied. On the other hand, this 
assumption allows us not only to simplify the proof of the existence of a weak solution and to consider 
more general functions / and b, but also to use the Moser-type iteration technique for improving the 
regularity of the constructed weak solution in the viscous case. 
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Once the existence of a solution is achieved, we show dissipativity of the (multivalued) dy- 
namical process associated to the system. This in particular entails existence of a (weak) trajectory 
attractor in the sense of Chepyzhov and Vishik [5] . 

Our subsequent results regard only the viscous case S > (and in particular can be applied to 
the Cahn-Hilliard model up to the modifications described above). Then, we can also prove that the 
trajectory attractor can be intended in the strong sense (i.e., w.r.t. the strong topology of the natural 
phase space) at least under slightly more restrictive conditions on /. The proof relies on an ad- hoc 
integration by parts formula and a variant of the so-called energy method (cf. [1] , see also [2l HH] and 
[TU1 [TT1 IT2"] for applications to trajectory attractors). If we furtherly restrict the class of admissible 
functions / (namely, asking n to be large enough) , then we can also prove that weak solutions become 
uniformly separated from for any time t > 0, so that the degenerate character of the system is 
actually lost for strictly positive times. This result, which is in our opinion the main achievement of 
this paper, is shown by means of a suitable version of the Moser iteration scheme which takes time 
regularization effects into account. As a further consequence of this "separation" property, we can 
also prove arbitrarily high regularity of weak solutions, as well as uniqueness, at least for t > 0. In 
turn, this permits to interpret the global attractor in the frame of the standard (single- valued) theory 
[U [55] , rather than in the trajectory sense. 

The plan of the paper is as follows. In the next Section [2] we will report our notation and 
hypotheses and the statement of our existence result. Its proof is divided into several steps and will 
be presented in Section [31 In Section [H we will show the existence of weak trajectory attractors. 
In Section [31 we will prove existence of a strong trajectory attractor in the viscous case S > by 
applying the so-called energy method. Finally, in Section [51 we will prove the strict positivity of u in 
the viscous case under more restrictive growth conditions on /. 



2 Existence result 

Let S7 be a smooth bounded domain of R d , d G {2,3}. Let T > a given final time and let Q :— 
SI x (0,T). Let H := J L 2 (S1), endowed with the standard scalar product (•, •) and norm || ■ ||. Let also 



V := 1 (SI) . We note by || • \\x the norm in the generic Banach space X. 
We make the following assumptions on data: 

b(r) = r s +/3r™, j3 > 0, r > 0; (2.1) 

< n < s < 10 if d = 3, < n < s if d = 2; (2.2) 

f(r) = ~, k>1, k>s + 1, r>0; (2.3) 

7 G VK 1 ^°°(R)n J L 1 (R); (2.4) 

geVnL°°(n), g:=\\g\\ L oo (a y, (2.5) 

5 > 0. (2.6) 



Remark 2.1. Our key assumption here is that the singular character of / has to dominate over 
the degeneracy of b at (cf. (|2.3[1 ). Actually, one could see with straighforward modifications in the 
proofs that in the case j3 > (i.e., if b(r) has a lower order degeneration r n at 0), then it would be 
enough to ask n > n + 1 rather than n > s + 1. We also point out that the requirement s < 10 in 
the three-dimensional case is motivated by the growth of b at oo (and not by its degeneration at 0). 
Consequently, in the application to the Cahn-Hilliard model (where solutions have to stay in between 
the two barriers r = ±1), s could in fact be arbitrary also in the case d = 3. 

We also define, whenever they make sense, the following functions: 

F(r) := -1_ + f(j) dr, F(r) := J 7 (r) dr, W(r) := F{r) + T(r). (2.7) 
Remark 2.2. According, e.g., to [B], a physically relevant expression for W' = f + 7 is given by 

W'{r) 77 + 4' where k G K )' ( 2 - 8 ) 
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Actually, this situation is not covered by our assumptions (I2.3[) - (|2.4[) . However, it is clear that, just 
with technical modifications in the proofs, one could replace f|2 . 3[) with something like 



ci^</'M<c 2 ^, (2.9) 

for all r > and some c\, C2 > 0. Assuming (|2.9p and properly choosing 7, it is clear that we can deal 
with the case (|2.8[) . Nonetheless, we will assume (|2.3I) in place of (|2.9I) in order to reduce technical 
complications in the proofs. 

Next, we introduce the energy functional associated to system (jl.ip - (|1.2l) : 

£{u):= f (^ + W{u)-gu). (2.10) 

This leads to defining the energy space, that will act as a phase space for our system: 

X := {u E V : u>0 a.c., u 1 '* £ L^O)}. (2.11) 

The space X is endowed with the natural (graph) metric 

d*(ui,tt 2 ) := IK -u 2 ||v+ |K~ K -«a~1Li(n)> (2.12) 

which is readily proved to be complete. Given u. € (0, 00), we also define 

%n ■— {u £ X : Un = fi}, (2.13) 

(•)n denoting here and below the spatial average over O. Actually, integrating (jl.ip in space one 
readily sees that the quantity un(t) is conserved in time for any solution u. 

Thanks to (|2.4p . (12.5P and to the above conservation property, a direct computation shows 
that, for some 0^,0^, > (also depending on g), there holds 

^(iNv + ||« 1_ 1Li (n) ) - ^ < £ ( u ) < c m(i + \Hv + h 1_ 1L 1( a))- ( 2 - 14 ) 

for all m E In particular, for a function it of assigned spatial mean /i, the finiteness of the energy 
£{u) corresponds exactly to the condition u E X^. 

The above notation is sufficient to define the class of weak solutions. 
Definition 2.3. A weak solution to problem (|l.l[) - (|1.4p is a couple (u,w), with 

u E L°°(0, T; V) fl L 2 (0, T; H 2 (fl)), 5 1/2 u e H 1 (0,T; H), (2.15) 

F(u) GL°°(0, TjL 1 ^)), (2.16) 

toG^'^Ti^p)), (2.17) 

fe 1/2 (u) Vw E L 2 (0, T; fl"), (2.18) 

6(ii)Vw6l*(g), m 4 E LnP(0,r ; (W" 1 '^)*^)), (2.19) 

/(«) E L 2 (0, T; L x (fi)) n L 5 / 3 (Q). (2.20) 

such that the following relations hold a.e. in (0,T): 

0) + (b(u)Vw, V(j>) = V^eW 1 '^^), (2.21) 

(w,iP) = 5(u u ii>) + (Vu,ViP) + (f(u)+ 1 (u)-g,iP) E V, (2.22) 

where (•, •) denote suitable duality pairings, and such that, in addition, 

u\t=o = uq, a.e. in ft. (2.23) 
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Remark 2.4. The above regularity framework actually refers to the three-dimensional case. If d = 2, 
(|2.17p - (|2.20[) could be improved a bit (we omit the details). 

We can now state our basic existence result, which can be considered as a variant of the theorem 
proved in [141 Sec. 8] (our assumptions are, indeed, slightly different). 

Theorem 2.5. Let us assume (|2.ip - (l2.6p and let, for some fi £ (0, oo), 

u £ X M . (2.24) 

Then, problem (jl.ll) - (|1.4|) admits at least one weak solution. 



3 Proof of Theorem [222 

We will give the proof in the case d — 3 and just point out some minor differences occurring in the 
two-dimensional case. 



3.1 Regularized problem 

First of all, we introduce a suitably approximated statement. Namely, given e £ (0, 1), we set 

6 e (r):=6((r a +e°) 1/a ), (3-1) 

{/(r) if r > e, 

/(e) + /'(e)(r -£) = -— + ^ TT (r-e) if r < e, 

where a > will be chosen later on, and it is intended that f e is defined for all rel. It is then worth 
noting that 

— rr if r > e, 
K(r) := <! ^ (3.3) 
if r < e. 



Moreover, setting 



we obtain that 



rK + l 



Fe(r) := -—^ + / Mr) dr, (3.4) 



if r > e, 



F e (r) := <( (K ^ 1 ! K (3.5) 

At this point, we can consider the approximate statement 

u e . t - div(6 £ (u £ )Vw e ) = 0, (3.6) 
w £ = Su e>t - Aw e + f e (u e ) + j(u e ) - g, (3.7) 

coupled with the initial conditions and the no-flux boundary conditions. Then, in analogy with 
[T31 Thms. 2 and 4] (see also [21 Thm. 2.1]), we have the following existence result for approximate 
solutions. 

Theorem 3.1. Let us assume (|2.1[) - (|2.6[) and let b e , f £ be specified by (|3.ip - (|3.2p . Let, in addition, 
u o,e S H 3 (£l) be defined as the (unique) solution to the elliptic problem 

u , £ - £ 2 Am , £ = iio, d n u Q . £ \ d Q = 0. (3.8) 

Then, there exists at least one couple (u e ,w e ) with 

u s £ L°°(0,T;V) nL 2 (0,T;# 3 (Q)), 5 1 ^ 2 u £ £ i? 1 (0, T: H), 5u e £ H 1 {Q,T;V), (3.9) 
w e £L 2 (0,T;V), (3.10) 
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satisfying (|X5 ]l -([57f )l a.e. in fi x (0,T), together with 



u E \t = o=u 0: e, a.e. in Q. (3.11) 
Moreover, in the case 6 > 0, tie couple (u e , w £ ) is unique. 

Sketch of proof of Theorem 13.11 We can proceed by following very closely the proofs of the 
quoted results [131 Thms. 2 and 4], [21 Thm. 2.1]. For this reason, we will just give some very brief 
highlights. Actually, the main difference here is due to the the growth of b £ at infinity. Nevertheless, 
one can of course truncate b £ near oo and replace it by some approximation b £V of at most linear 
growth, such that b £l , suitably tends to b e as v \ 0. Then, existence for v > is proved similarly as 
in [131 13] and it remains to prove suitable a-priori estimates uniform w.r.t. v. However, for the sake 
of simplicity we will omit the ^-approximation and rather perform formal estimates on the e-solution 
in order to show that it fulfills regularity properties ()3.9|) - (|3. 10|) . 

Thus, we can firstly perform the energy and entropy estimate, as below. Using that b £ > e as / 2 
and the Lipschitz continuity of f £ , we then obtain 

\\ u s\\L^{0.T;V)nL- 2 (0,T:H 2 (n)) + <5 1/2 ||We||H 1 (0,T;H) < Ce, (3-12) 

\\bl /2 Vw £ \\ L 2(0,T;H) + \\w £ \\l 2 (o,T;V) < c £ . (3.13) 

Next, we test (|3~7)) by A 2 u £ . Using (|3~8|) . f[37T2]) - ([37TB]) . ([2~5j) and the global Lipschitz continuity of f £ 
and 7, it is then not difficult to obtain 

IM|£2( 0) T;ff3(fi)) + S 1/2 \\u £ \\ L °°(0,T;H 2 (n)) < Ce- (3-14) 

In case 6 > 0, we can also test (|3.7p by — Au e . t , that yields 

tf||Ue||ffi(0,T;V) < Ce- (3-15) 

This gives all desired regularity properties. Finally, to prove uniqueness in the case S > 0, we can test 
the difference of (|3.7[) by the difference of the u £j t and the difference of (|3.6[) by the difference of w £: t- 
The details are left to the reader. I 



3.2 A priori estimates 

We now aim to obtain a number of a priori bounds, uniform in s, with the purpose of removing the 
approximation. 

Energy estimate. We test (|3.6|) by w e , (13.7[) by u e t, and sum the results. We then obtain 



where 



-£ £ (u £ ) + S\\u £ , t \\ 2 + / b £ (u £ )\Vw £ \ 2 = 0, (3.16) 
Jn 

£ £ {u) := / (JM + Fe{u) + T (u) gu) . (3.17) 

Remark 3.2. We point out that, even at the approximate level, this formal estimate is not completely 
justified in the case 8 — 0. Actually, w £ is (only) in L 2 (0,T;U), while (I3.6[) is not an equation in 
L 2 (0, T; V') since b £ grows fast at infinity. However it is clear that, performing a truncation of b £ and 
then passing to the limit, the estimate could be justified. 



Then, we integrate (I3.16[) in time and notice that, by (|3.5p . 

F e {u , e )< [ FM< f F(«o) + (/(uo, e ),uo, e -uo) < / F(u ), (3.18) 
: Jn Jn Jn 

the latter inequality following from (13.81) and Green's formula. Thus, owing to (|2.24[) . (|3 . gives 

IKIU°°(o,T;V) + ^ 1/2 ||u £ ,t|U 2 (o,T;/f) < c, (3.19) 

\\F e (u e )\\ L oo (0!T . L r m <c, (3.20) 

|j6i/> £ )V W£ || L2(0 , T;H ) <c. (3.21) 
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Here and below, c denotes a positive constant, independent of e and of time, whose value may vary 
even inside a single row. We will use the letter a to denote constants used in estimates from below. 

Entropy estimate. Let us define, for r G (0, oo), the entropy M, by setting 

m ( r ) := / ITT' M W := / m ( r ) dr - ( 3 - 22 ) 



b(r) 

Clearly, M is a convex function that grows at most like r for r ~ oo. Moreover, let us introduce its 
approximate version by taking, for r G K, 

m e (r) := jf ^y, M £ (r) := ^ m £ (r) dr. (3.23) 

Clearly, M £ is a convex function such that M £ < M a.e. in (0, oo). Moreover, m £ , M £ tend to m, M, 
respectively, uniformly on compact sets of (0, oo). 

Then, we can test (|3.6I) by m s (u e ), (|3.7|) by — Au £ , and sum the results. We deduce 

±([ M £ (u £ ) + 6 -\\Vu £ \\ 2 )+\\Au £ \\ 2 + f (f>(u £ )+y(u £ ))\Vu £ \ 2 + (g,Au £ )=0, (3.24) 

and we have to control some terms. Firstly, by (|2.4p and Holder's inequality, we have 

7 > £ )|Vu £ | 2 + (g,Au e ) <c(l + ||V Ue || 2 ) + i||A Ue || 2 . (3.25) 

We now observe that, for all e e (0, 1), it is M £ (uo t£ ) < M(u ) (to prove this, proceed as in (|3.18t ). 
Moreover, we have that M(u ) G L 1 (fi). Actually, (|2.24l) entails G L j (0) and M(r) grows no 

faster than r 2 ~ s in the neighbourhood of 0, which is good since we assumed k > s + 1 (cf. (|2.3p ). 
Thus, integrating p.25p in time we arrive at 

\\ u e\\L 2 (0,T-H 2 (n)) < c, (3.26) 



T 

' f £ (u £ )\Vu £ \ 2 <c. (3.27) 
o Jn 



Control of the nonlinear terms. Here, our aim is to derive e-uniform bounds in order to pass to 
the limit in the terms b £ {u £ ) and f £ (u £ ). 

First of all, by (|3.19[) . (|3.26p and interpolation, we obtain, if d = 3, 

||We|U 10 (fix(0,T)) < C, (3.28) 

whereas for d = 2 we have instead 

||«e|Uf(nx(o,T)) < cp Vp G [1, oo). (3.29) 

Let us now set 

Ql := {(x,t) G Q : u £ (x,t) < s}, (3.30) 

Ql := {(x,t) G Q : e < u s {x,t) < l}, (3.31) 

Ql := {(x,t) G Q : u £ (x,t) > 1} (3.32) 

and notice that, by (|3.5[) and (|3.20p . 

\ttl(t)\ < cs"" 1 for a.e. t G (0,T). (3.33) 

Here and below, fi*(t), i G {1,2,3}, is the section of Q\ at the generic time t G (0,T). Also, it is 
obvious that 

IIA(«e)||i-(03) <c. (3.34) 
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Next, we remark (cf. (JH3J) that 

1 f —p, — if r > e, 

Mr)-1 « if r < e . (3 ' 35) 

V. £-as/2 

Hence, using p. 51) . (|3.20|) . (|3.33|l and the condition k > s + 1 (cf. (|2.3|l ). taking a small enough in the 
definition (13.11) of 6 e we obtain 

lk _1 ML~(0,T;Li(n)) ( 3 - 36 ) 

Setting now 

z e (x, t) := max {u e (x, t), e}, ti e (x, t) := z e (a;, i) -5- , (3.37) 
for a.e. £ £ (0,i) we have 

ll«e(t)|l!a (fi) <c(||« e Wf + ||V« e (t)|| 2 ) 

<c/ A + c/ 4t+ C / ^ (3-38) 

(here and below, we compute the exponents referring to the case d = 3; for d = 2, they can be 
improved, of course). Hence, using (|3.20[) . (|3. 241) and (|3.33l) . and integrating in time, it is not difficult 
to arrive at 

II^IU 2 (0 : T;LS(f2)) < C, (3.39) 

whence, recalling (I3.35[) and possibly choosing a smaller a, 

\\ b e\u s )\\ L1{QtT . L3m) <c. (3.40) 
Then, using (I3.2ip and either (I3.36|) or (|3.40[) . we also have 

||Vw £ || L 2 (0!T;i i (o)) + ||Vu> e ||£i( 0)T; £3/2( n )) < c, (3.41) 

whence, by interpolation, 

l|Vu? e ||£V4(Q) + l|Viu £ ]| I ,5/ 3 ( 0)T . z , 1 5/i4(n)) < c. (3.42) 

We are now ready to give an estimate of the term f e (u e ). To do this, we first test (|3.7|) by u e — fx, to 
obtain 

A* / \fe(ue)\= / ( ~ <5u M + Au £ - j(u £ ) + g)(u e - fi) + / w e (u e -fi)- / f s (u e )u E (3.43) 
Jsi Jn Jn Jn 

and the terms on the right hand side are treated as follows: 

( - 5u Sjt + Au £ - j(u e ) + g)(fi - u s ) < c\\ - 5u £:t + Au e ~ j(u s ) + g\\\\fj. ~ u e \\ 

< c\\ - Su e , t + Au e - 7(u e ) + g\\ =: 771, (3.44) 
thanks to ([53U)> . where ||j7i||z,2( ,t) < c by ([231), and (E^SJ). Next, 

;(u E -fx)= (w e - (w e )a)(fx-u E ) < c||Vto e || =: (3A5) 
Jn 

with ||t72||l 2 (o,t) < c by (|3.41[) . Finally, due to p. 21) it is clear that 

fe(Ue)u s < § / |/ e (n e )|+C M . (3.46) 

Squaring (|3.43p . integrating in time, and using (|3.44p - (|3.46l) . we easily arrive at 

||/ e (u e )||z,2(o,T;Z,i(n)) < c (3.47) 



We 

u 
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Setting now <j> e := — |/ e (r)| 2 / 3 , we test (|3.7p by 4> e (u e ) — ((j)e{u s ))a, and integrate in space and 
time. This gives 



o Jn 

T 



fe{u e )<l>e{us) + / (Su e , t + "{(u s ) - g,<f> e (u e ) - ((f> e (u E ))n) 
Jo 

(w s ) Q ^ e (u s )) + f (i<k(ue)h [ A(«e))« (3.48) 



'o Jo 
Let us first notice that the first term on the left hand side gives 

rT 



f e (u e )<t> e (u e ) = \\f e {u e )\\% {Q y (3-49) 



io Jn 

Next, by Holder's inequality, we obtain 

rT 



(Su E:t + 7(u e ) - g, <f) s {u e ) - (<t) e (u e ))n) 
< \\Su Stt + j(u s ) + g\\ 2 L 2 (Q) + c\\<j> e {u e )\\ 2 L2[Q) < <j\\fe{u £ )\f/? /3{Q) + c<r- (3.50) 

where a is a small constant to be chosen below and c a > depends on a. Notice that (|3.19[) . (|2.4[) 
and (|2.5[) have been used here. 

It then remains to control the terms on the right hand side of (|3.48|) . As far as the first one 
is concerned, recalling (|3.34p . we notice that, a.e. in (0,T), 

(w £ - (w E )n,<l> e (u e )) < \\cj) e {u e )\\ L 5/2( Q) \\w e - (w e ) n \\ L 5/3 {Q) 

< ^\\fe(Ue)tJs/3 (Q) + CaWVw^l,^. (3.51) 

Finally, let us estimate the latter term in Q3.48p . Actually, it is clear that 



{4>e{u e ))u / fe(u e ) < 4 f £ {u e ) || L 2 (0iT:L i (n )) + C < C, (3.52) 

Jn 

the latter inequality following from (|3.47p . Collecting now (13.481) - (I3.52p and recalling (|3.34p . we finally 
arrive at 

||/e(«s)|| i «/.(Q)<C (3.53) 

3.3 Passage to the limit 

We will just consider, for brevity, the case d = 3. For simplicity of notation, let us set £ £ := 

1 /2 

b e (u £ )Vw £ and let || • || p denote the norm in the space L P (Q). Then, by (|3.2ip . 

HCella < c (3.54) 
Moreover, being s < 10, (|3.28|) guarantees that 

||& e (u E )Vw e || = ||6y 2 (M e )C e || < c for some q > 1. (3.55) 
Hence, by comparison in (|3.6[) . 

||«e,tlU«(0,T;W-i.<i(n)) < c for some q > 1 (3.56) 

(of course, if d > we have much more, but we want to deal with the most general case here). 
Consequently, using (|3.19p . (|3.28[) . the Aubin-Lions lemma, and Lebesgue's theorem, 

u e ->• u strongly in L q {Q x (0,T)) Vq G [1, 10) (3.57) 
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(here and below, all convergence relations arc to be intended up to the extraction of non-relabelled 
subsequences). Let us now notice that 

b e — > b uniformly on compact subsets of R, (3.58) 

where b denotes the even extension of b to R. Then, by p.57[) . s < 10 in (|2.ip . and Lebesgue's theorem 
again, we obtain 

b e (u £ ) -> b{u) strongly in L q (n x (0,T)) for some q > 1. (3.59) 

Analogously, by (|3.53D . 

f e (u e ) -> f(u) strongly in L q (Q x (0,T)) for all q e [1,5/3), (3.60) 

whence, in particular, the limit u is a.e. nonnegative and we can replace b with b in (|3.59p . 

Our next aim is to pass to the limit in the product b e {u E )Vw e . To do this, we first notice 
that, by (|3~57)l and s < 10, 

b £ (u E ) 1/5 ->• b{u) 1/5 strongly in L q {Q) for some q > 5. (3.61) 

Thus, using also (|3.42[) we arrive at 

b e (u e ) 1/b Vw e -> 6(u) 1/5 Vw weakly in L q (Q) for some g > 1. (3.62) 

Next, interpolating between (13.361) and (|3.40p . we get 

||6 s (« e ) _1 || B/3 < c. (3.63) 

Using (|3.63p and (|3.54p . we then obtain 

||& e ( U£ ) 1/5 V We || 25/17 < ||C £ || 2 ||6 £ (^)- 3/10 || 50/9 < c. (3.64) 

Hence, by (pTM)) , (|3~HT1) and 

i £ (« E ) 2/5 Vw £ = 6 B (u e ) 1/B (6 e (us) 1/B Vu; e ) -> 6(u) 2/5 Vw weakly in L 25 / 22 (Q). (3.65) 

Then, writing 6 e (w e ) 2//5 Vu; e as Ce^e( M e) -1 ^ 10 and using once more (|3.54l) and fl3.63[) . it is clear that 
the exponent 25/22 in (|3.65[) can be improved. Thus, iterating the above procedure and using once 
more (|3.54p . it is not difficult to arrive at 

Ce^b{u) 1/2 Vw weakly in L 2 (Q), (3.66) 

whence, using once more (|3.59p (and computing explicitly the exponent), we finally obtain 

b e {u e )Vw e ^b{u)Vw weakly in LW^(Q). (3.67) 

Thus, we can take the limit of all terms in (|3.6p - (|3.7D and get back (|2.2ip ~ (|2.22p . where the properties 
required to test functions depend of course on the regularity conditions proved above. Indeed, (|2.15l) - 
(|2.20p follow as a direct byproduct of the procedure. The proof of Theorem 12.51 is concluded. 

4 Weak trajectory attractors 

In this section, we construct the so-called weak trajectory attractor for problem fj 1 . 1 11 - (11 .2[) . To avoid 
technicalities, we will limit ourselves to deal with the (more degenerate) case j3 = 0, i.e., b(u) = u s , 
for s e [0,10) (cf. flUD ), 

We start by proving a dissipativity result holding for the weak solutions constructed in the 
proof of Theorem 12.51 
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Theorem 4.1. Let the assumptions of Theorem 12.51 hold, with (3 — 0, and let 

E a :=£(u ), (4.1) 

that is finite thanks to (12.241) . Then, there exist a weak solution (u,w) and a monotone function 
Q : [0, oo ) — > [0, oo) such that 

£(u(t))+ [ (5\\u t \\ 2 + \\b 1/2 (u)Vw\\ 2 ) <Q(E ) Vi>0. (4.2) 



More precisely, there exists a set Bq bounded with respect to the metric (I2.12|) . such that, for any 
dx -bounded set B C X^, there exist a time Tb > such that for any initial datum u e B there 
exists at least one weak solution u starting from ii and such that u{t) £ Bq for all t >Tb- 

Remark 4.2. Notice that the above result does not claim that dissipativity holds in the whole class 
of weak solutions, but just that from any admissible initial datum there starts (at least) one weak 
solution in the dissipative class (cf. Remark 14.71 for further considerations). 

Proof. We integrate (|3.16l) between and an arbitrary t > 0. This gives the e-equivalent of (|4.2p . 
Then, we take the liminf with respect to e \ and use estimates (|3.19[) . p.57[) . (|3.66[) . the fact that 
F £ converges to F uniformly on compact sets of (0, oo), Fatou's Lemma, and the lower semicontinuity 
of norms with respect to weak or weak star convergences (of course, we will get an energy inequality, 
and not necessarily an equality, in this way). Relation (|4.2I) is proved. 

To show dissipativity, we start considering the case k > s + 1. Then, it is convenient to rewrite 
the energy inequality in the differential form: 



—8{u) + 5\\u t \\ z + / b{u)\Vw\ z <Q, (4.3) 

for a.e. t > 0. We then notice that, setting z := u^ 1 , we have, from (|2.14l) . 

viW^hw + llVuf) - c < £{u) < a~ l H|z||2„ (n) + ||Vu|| 2 + l), (4.4) 

where a G (0, 1) depends in particular on fi = uq. 

Now, similarly with (|3.43[) . we test (|1.2I) by u — (i, obtaining 



||V«|r + /i||2||2» ( n)< / (Su t -j{u)+g){u-fi) + \\z K L \\ L i m + / w(u-n) 
Jn Jn 

< i||V U || 2 + c(l + S 2 \\u t \\ 2 ) + ^\\z\\ K LK(n) +c^ + c\\w - wn\\ 2 L3/2m 

< I|| Vu || 2 + c(l + 5 2 \\u t \\ 2 ) + |||z||2» ( o) + + c\\Vw\\l Ha) . (4.5) 
The last term can be controlled this way: 

\\Vw\\ 2 L1(n) <c\\z s \\ LHn) [ b(u)\Vw\ 2 . (4.6) 
Jn 

Next, being k > s+1, using the first inequality in (|4.4[) and Jensen's inequality, we obtain, for suitable 
positive constants a, 

(£(u))^ > a\\z s \\ Lim > a{z n ) s > a(u n )- s = a^ s . (4.7) 
From (|4.4p - (l4.7p . we then obtain, for a suitable c* > 0, 

£{u) <c*(l + 5 2 \\u t f + (£(u))^ J b(u)\Vw\ 2 y (4.8) 

Thus, assuming £ > 1, which is of course not restrictive, we can divide by £ S /( K_1 ) to obtain 

<c,(l + 5 2 || Ut || 2 + f b(u)\\7w\ 2 y (4.9) 
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Taking the (1 — e)-power for e £ (0, 1), summing to (14.31) . and applying Young's inequality, we then 
have, for some a > 0, 

^£(u) + a((£(u)) ( "" 1 »-i 1 " e> + s\\u t \\ 2 + J b(u)\\7w\ 2 ) < c e , (4.10) 

whence the thesis follows by integrating in time and applying the comparison principle for ODE's. 
More precisely, we also have a quantitative decay estimate for the energy. 

In the case k = s + 1, we can say a little bit less, but dissipativity still holds. Actually, we 
can repeat the above procedure up to (|4.8p . Then, we notice that (|4.3p implies in particular that 



(5\\u t \\ 2 + J b(u)\Vw\ 2 ^ < E < oo. (4.11) 
Consequently, there exists at least one time T* such that 



T* G [0,2c*E ], and 5\\u t (T*)\\ 2 + j b(u(T*))\\7w(T*)\ 2 <—. (4.12) 



Substituting in (|4.8[) . we then have 

£(u(T,))<2c,+S=:C*. (4.13) 

Then, we obtain £(u(t)) < C* for any t > T* simply observing that, by (I4.3[) . £ is nonincreasing. I 

As a next step (see [5] for more details), we need to rewrite the dissipative estimate in such a way 
that, on the one hand, we will be able to control all the norms which are necessary to pass to the 
weak limit in the space of solutions of the problem considered (and verify that the limit function is 
again a solution) and, on the other hand, be sure that the corresponding trajectory phase space will 
be translation invariant. 

The following lemma improves the dissipative estimate (|4.2p by adding the terms controlled 
by the entropy estimate. 

Lemma 4.3. Let the assumptions of Theorem 14.11 hold. Then, there exists a solution (u, w) of 
problem (|l.ll) - (|1.2p which satisfies the following estimate: 

£{u{T)) + [ T+1 5\\u t (t)\\ 2 + ||6 1 / 2 ( u (t))V«;(i)|| 2 + \\Au(t)\\ 2 + (/'(«(t))V«(t), V«(t)) dt 
Jt 

< Q(£(u(0)))e- aT + C*, VT>0. (4.14) 

where the positive constants a and C» and the monotone function Q are independent of t and of the 
concrete choice of the solution u. 

Proof. As a consequence of Theorem l4.11 it is clear that (|4.14p holds, for suitable Q and a, without 
the last two terms in the integral on the left hand side. To control these terms, it is sufficient, in the 
case k > s + 1, to sum (I3.24p to (|4.10p in the preceding proof. In the case K = s + 1, we know that 
there exists T* = T»(Eo) such that the energy is smaller than some constant C* independent of the 
initial data for any T > T* (cf. (j4"T3l ). Then, it is sufficient to integrate ([3T24]) over (T, T + 1) for 
T > T* to get 

T+1 ||Au(t)|| 2 + (/'(u(t))Vu(t), Vu(t)) dt 

T 

<c + J^M(u(T))+ 5 -\\Vu(T)\\ 2 < Q(£(T)) < Q(C*), (4.15) 

as desired. Actually, it is clear that the energy controls from above the terms M(u) and ||Vw||. Being 
pedantic, all these estimates should be done on the level of approximations u e with passing to the 
limit after that; we directly performed the estimate on u just for brevity. I 
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Remark 4.4. Note also that all the norms involved into our definition of a weak solution (see Theorem 
I2.5[) are under control if we assume that the weak solution satisfies (|4.14p . This fact, which can be 
verified exactly as in the proof of Theorem [531 is crucial in order to be able to pass to the weak limit 
on the space of weak solutions and establish that the absorbing set for the trajectory semigroup is 
indeed closed, see below. 

We are now able to define the trajectory phase space and trajectory dynamical system associated with 
problem (fTT>([TT2)l . 

Definition 4.5. Let K.+ C L°°(M. + ,£) be the set of all solutions u of problem 1|1.1[1 - (|1.2|) belonging 
to the class (|2.15p - (|2.20[) which satisfy the following analogue of (|4.14|) : 

£(u(T)) + [ T+1 S\\u t (t)\\ 2 + ||&Va( tt ( f ))v«,(i)|| a + \\Au(t)\\ 2 + (f(u(t))Vu(t), Vu(t)) dt 
Jt 

<C u e~ aT + C*, VT>0. (4.16) 
for some constant C u depending on the solution u. Then, the shift semigroup T(h), h > 0, acts on 

T(h) : K.+ -> JC+, (T(h)u)(t) := u(t + h). (4.17) 

We will refer below to 1C+ and T(h) : /C+ — > /C+ as a trajectory phase space and trajectory dynamical 
system associated with problem (|l.l|) - l|1.2|l respectively. 

Furthermore, in order to be able to introduce the attractor of the trajectory dynamical system, 
we need to specify the topology on /C+ as well as the class of bounded sets. 

Definition 4.6. We endow the set K.+ with the topology induced by the embedding K.+ C 0™ eak ~ 
[L£ c (R+,H 1 (n)) n Lf oc (R + ,H 1 {fl))] w * , where w* stands for the weak-star topology and will refer to 
it as a weak topology on the trajectory phase space JC+ . 

A set B C /C+ will be called bounded if inequality (|4.16[) holds uniformly with respect to all 
u € B, i.e., if 

C B ■= sup C u < oo. (4-18) 

uEB 

Remark 4.7. As usual (see [9] for the details), under the general assumptions of Theorem 14.11 we 
know neither the fact that any weak solution of problem (ll.ll) - (|1.2j) satisfies the energy inequality 
(|4.16p nor that the constant C u in (|4.16p can be expressed in terms of £(u(0)). Actually, it may be 
possible to construct a solution u which satisfies (14. 14[) for the initial moment T — only and be unable 
to verify its analogue for other initial times. By this reason, attempting to replace (|4.16p by (|4.14p 
in the definition of the trajectory phase space /C+, we lose the translation invariance T(h)K,+ C /C+ 
which is crucial for the attractors theory. However, as we will see below, under the more restrictive 
assumptions of Theorem l5.11 the answer to both the questions posed above is positive. So, in that case, 
every reasonably defined weak solution satisfies (|4.14[) and the boundedness condition is equivalent to 
the boundedness of w(0) in the energy space. 

Finally, we are now able to introduce the trajectory attractor for problem (|1.1|) - (|1.2|) . 

Definition 4.8. set A tr C /C+ is a (weak) trajectory attractor for problem p.l[) - (|1.2[) if the following 
conditions are satished: 

1) A tr is compact in e^ eafe ; 

2) It is strictly invariant with respect to the trajectory semigroup: T(h)A tr = A tr ; 

3) It attracts the images of all bounded sets of /C+ as time tends to inhnity, i.e., for every 
bounded subset B C K.+ and every neighborhood 0(A tr ) (in the topology of 0T eafe / ), there exists a 
time T = T(B, O) such that 

T(h)B C 0{A tr ), 

for allh>T. 

Next, we can state the existence result for the above introduced object. 
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Theorem 4.9. Let the assumptions of Theorem 14.11 hold. Then, problem (jl.ip ~ p.2p possesses a 
trajectory attractor A tr in the sense of the above definition. Moreover, this attractor is generated by 
all complete (i.e., defined for all t g M.) and bounded trajectories for that system. Namely, we have 

A tr :=K\ t > , (4.19) 

where JC C L°°(R, £) is the set of all solutions of (fl~T j) - (fT2")) which satisfy 

E(u(T)) + [ T+1 S\\u t (t)\\ 2 + H^K^V^f + ||A U (<)|| 2 + (f(u(t))Vu(t), Vu(f)) dt < C„ 



JT 

for all T G R and some C* > 0. 

Proof. As usual (see in order to show the attractor existence, we only need to verify the 
existence of a compact and bounded absorbing set for the trajectory dynamical system T(h) : JC+ — > 
JC+ (the continuity of the semigroup in the 0™ eafe -topology is obvious since it is just a translation 
semigroup). Note that, due to (|4. 16[> . the set B C K.+ of solutions u satisfying 

E(u(T)) + J™ S\\u t (t)f + \\b^ 2 (u(t))^w(t)f 

+ || Au{t)\\ 2 + (/'(u(t))Vu(t), Vu(t)) dt < 2C*, (4.20) 

for all T > will be an absorbing set for the semigroup T(h) acting on /C+. Obviously, this set is 

bounded (in the sense of the Definition I4.6|> . Thus, we only need to verify that it is compact in the 
Qweak topology 

Indeed, let {u n } C B be a sequence of solutions. Then, due to estimate (|4.20l) . this sequence 
is precompact in 0™ eafc ; so, without loss of generality, we may assume that u n — )• u € 0™ eak in the 
topology of 0™ eafc and we only need to verify that the limit function u solves (|1.1|) - (|1.2|) and satisfies 
(|^2H)l as well. 

The proof of this fact repeats almost word by word the proof of the existence Theorem 12.51 
and is even a bit simpler since we do not need to consider the regular approximations to / and b 
(note that the uniform estimate (|4.20p allows us to control uniformly all of the norms involved into 
(|2.15p - (|2.20p ). By this reason, we leave the rigorous proof to the reader. 

Thus, all of the assumptions of the abstract attractor existence theorem are verified and the 
theorem is proved. I 

5 Energy equalities and strong attraction 

In the viscous case S > and under slightly more restrictive assumptions on the growth of /, we 
can prove that At? is in fact a strong trajectory attractor (i.e., it attracts with respect of the strong 
topology of X). This is the object of our next result: 



Theorem 5.1. Let assumptions (12.1l) - (|2.6p hold and let, additionally, 5 > and {3 = 0. In addition, 
let 

3s 

k > — — h 1 if d = 3, and k > s + 1 if d = 2. (5.1) 
Then, any weak solution of problem (|l.ip - (ll.2p satishes the additional regularity properties 

w e L 2 (0, T; H), f(u) G L 2 (0, T; H). (5.2) 
Moreover, a.e. in (0, oo), there holds the following energy equality: 

^-£(u) + 6\\u t \\ 2 + [ b(u)\\7w\ 2 = 0, (5.3) 
dt Jn 

as well as the following entropy equality (compare with (|3.24p ): 

1( j M(u)+ 5 -\\Vu\\ 2 )+\\Au\\ 2 + ( (f(u) +1 i(u))\\7u\ 2 + (g,Au)=0. (5.4) 
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Proof. Let us start proving (|5.3p and first deal with the 3D-case. The key step is given by the 
following integration by parts formula. 

Lemma 5.2. Let b G L p (fl) for some p > 1, with b > a.e. in f2. Let also b^ 1 E L q (Vl) for some 
q > 3/2 if d = 3 (respectively, for some q > 1 if d = 2). Let <fi E H and let w be the (unique) solution 
to the degenerate elliptic problem 

- div(fcVw) + w = <f>, in fl, (bVw) • n = 0, 021 T. (5.5) 

Then, b\\7w\ 2 G L x (fi) and 

( - div(6Vw),w) = / 6|Vw| 2 . (5.6) 

J Q 

Remark 5.3. As it will be clear from the proof, the regularity of w is sufficient to state (|5.5p in that 
"strong" form. In particular, since we have 

bVw G L^+t(Q), div(6Vw)eff, (5.7) 
a suitable trace theorem (cf., e.g., [_ Thm. 2.7.6]) permits to interpret the boundary condition in 

_ p+l _2p_ 

(|5.5p in the sense of trace operators as a relation in the space W 2 p 'P+ 1 (r). 

Proof of Lemma 15.21 Again, we prove the theorem for d — 3 and just point out some minor 
differences occurring for d = 2. Let A be the Laplace operator with 0-Neumann boundary conditions, 
namely, 

A-.V-tV', (Av,z) := [ Vu • Vz. (5.8) 

Then, one can see A + Id as a strictly positive unbounded operator on and consider fractional 
powers of it. For s > 0, we let b £ :— max{6, e}. We now consider the approximate problem 

eA 3 w e - div(b £ Vw £ ) + w e = (j), (6 £ Vui £ ) ■ n = on T. (5.9) 

Then, testing (|5.9p by w e we obtain 

eP 3/ Vf+ / b s \Vw £ \ 2 + \\w e \\ 2 = (^w £ ). (5.10) 

Thus, for all e > 0, we have that w £ G L>((A + Id) 3 / 2 ) C H 3 (n). 

From (|5.10[) . we obtain that w £ is bounded, independently of e, in ff. In addition, we have 

IIV^I^^ < ||6^ 2 V W£ ||||6 e - 1/2 || L2 , (0) < c (5.11) 

and, being q > 3/2, it follows 2q/(q + l) > 6/5 (respectively, if d = 2, from q > 1 we have 2q/(q+ 1) > 
1). Thus, by standard compact embedding results, we have, up to a (nonrelabelled) subsequence of 
e\0, 

w £ — > w weakly in W (fi) and strongly in i/. (5-12) 
Moreover, being p > 1, we can write 

II^V^H < ||^/ 2 V W£ ||||6y 2 || i2p( o) < c. (5.13) 

LP+ 1 (ft) 

Thus, 6 e Vw E is bounded in L 2p ^ p+1) (f2) C D((4 + Id)" 1 ) and, consequently, 

|| - div(6 e Vw £ )|| Z5((j4+w) _ 3/2) < c (5.14) 



and, proceeding similarly with Subsection 13.31 we can also prove that b £ \7w £ tends to b\7w weakly 

2p 

in Lp+i^). Moreover, (I5.14p tells us that, for any e > 0, equation (|5.9p makes sense at least as a 
relation in D((yl + Id)~ 3 / 2 ) (in particular, the obtained regularity w £ G D((A + Id) 3 ^ 2 ) justifies having 
used w £ as a test function in (|5.9[l ). 



15 



Now, the obvious fact that b e —> b strongly in L p (£l), the first of (|5.12p . and Ioffc's theorem 
(see, e.g., [T7]) give 

/ b\Vw\ 2 < liminf / b e \Ww e \ 2 . (5.15) 
Thus, using (|5.10l) and the second of (|5.12[) , we can go on as follows: 



/ 

■hi 



b\Vw\ 2 < lim(0- iu e ,iu e ) + liminf ( - e\\A 3/2 w £ \\ 2 ) 



< 



w, w) = ( — div(bVw), w) , (5.16) 



where (|5.5p has been used to deduce the last equality. Thus, to complete the proof, we have to show 
the inequality converse to (jSTTCl) . Now, bVw £ L 2v ^ v+V > (fi) thanks to (the liminf of) ([535]) . Thus, 
also - div(6Vw) € £>((^4 + Id)~ 3/2 ) so that we can test ([53)1 by w e £ D((A + Id) 3/2 ) and rigorously 
integrate by parts to obtain 



<wJ 2 



(-div(bVw),w e ) = { bVw-Vw e <~ j b\Vw\ 2 + ~ f b\\7i 
< 1 /" 6|Vu;| 2 + i / 6 £ |V W£ | 2 

1/^2 J £"2 

-5 / fi & l Vw l 2 -fll^ 3/ Vl! 2 + i(0-w £ ,^), (5.17) 

where the fact that b < b e almost everywhere and the equality (|5.10l) have also been used. Passing to 
the limit and using the second (|5.12l) and that — div(6Vw;) £ H (as it follows by comparison in (|5.5I) ). 
we then obtain 

If 1 If 1 

( - div(6Vw),w) < - / b\Vw\ 2 + -{<j> - w,w) < - / b\Vw\ 2 + - ( - div(bVw), w) . (5.18) 

Namely, we obtained the inequality converse to (|5.16l) . whence the thesis. I 

We now proceed with the proof of Theorem 1 5. II and, precisely, of equality (15.31) under the assumption 
(|5.1[) . To do this, we first prove (I5.2[) and, with this purpose, we set z := it" 1 and observe that 
equation (11.21) can be rewritten as 

Sz t + z 2 Az~ 1 + z K+2 = -~z 2 w + z 2 <f>, where <\> := 7(1*) - g (5.19) 

and we notice that 

IMU~(nx<o,T)) < C, (5.20) 

thanks to (|2.4[) - (|2.5p . Here and below, C denotes a constant possibly depending on the "energy" of 
the initial data (cf. (|2.24p ) and on the choice of T, while c is an absolute constant (i.e., it does not 
depend on the initial data or on T). By (|2.18p . we also know that 

\\u s/2 Vw\\ LHnx(0X)) <C. (5.21) 

Next, by the first of (|2.20p . a comparison in (ll.2p gives also 

HHU»(o,T;£i(n)) < O. (5.22) 
At this point, we note that, for d = 3, thanks to (|2.16l) and (|5.ip . 

||&- 1/a (u)IU~ ( o lT ;za(n)) < C. (5.23) 
Hence, combining (|5.23p with (|5.22p and (|2.18p . we readily arrive at 

IIHU 2 (Q) < c \\ w \\L^(0,T-W^e/5 {n)) < C, (5.24) 

whence (|5.2p follows simply by comparing terms in (|1.2p and taking advantage of (|2.15p ~ (|2.20p . In 
the case d — 2, we only have the L°°(0, T; iJ)-norm in (I5.23p . but (|5.2p still follows by using the 
continuous embedding W 1 ' 1 ^) C i? in the analogue of (|5.24p . 
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We now proceed with the proof of (|5.3[) . Summing together (jl.ip and (|1.2p . we have 

- div(6(w)Vw) + w = (f> := (S - l)u t - Au + /(«) + -f(u) - g (5.25) 

and it is clear from (HO)) - (|2~5j) . (l2~l~5)) and (|5~2l) that, for any T > 0, it is ||^|U2( 0)T;iJ -) < C T . 

Moreover, from ([2^2]) and (the limit of) (l3~28)l (or (ET2910 . we have that, a.e. in (0,T), b(u) G 
L P (Q) for a suitable p > 1 (e.g., if d = 3, we can take p = 10/s). Finally, thanks to (|5.2p . it is clear 
that (if d = 3, the case d = 2 being analogous), a.e. in (0,T), 

||& -1 (t*)IU<(n) < W^h^n) < |U. ( n) + c < C, (5.26) 

for a suitable g > 3/2. Hence, b = b(u) and 4> defined in (|5.25[) satisfy, a.e. in (0, T), the assumptions 
of Lemma T5. 2 1 Consequently, as we test (jl.ll) by w and (|1.2j) by Uj, we obtain, thanks to (|5.6[) . the 
energy equality (I5.3f) . as desired. 

Finally, we come to the proof of (I5.4j) . From (|2.20[) and (|5.ip . we see that at least m(u) G 
L 2 (Q), so since due to (12.151) u t G L 2 (Q) as well, we have 

^(Af(«),l) = («t,m(u)) = (div(6(u)Vw),m(u)), (5.27) 

where the right-hand side is understood as a scalar product in L 2 (Q). Thus, we need to verify that 

(div(b(u)Vw),m(uj) = -(b(u)Vw,Vm(u)) = -(Vw,Vm) = (w,Au), (5.28) 

almost everywhere in time. To this end, we note that, according to (|2.15l) and the embedding H 2 (il) C 
C(H), we have u(t),b(u(t)) G C(fi) for almost all t. Then, keeping in mind that 6 1 / 2 (u)Vw G L 2 (Q) 
by (|2TT5j) . we conclude that 6(u(t))Vtfl(t) G i 2 (fi) for almost all <. 

Thus, we only need to check that m(u(t)) G i? 1 (fi) for almost all The fact that this function 
belongs to L 2 (il) is already verified, so we need that 

Vm(u(t)) = 6" 1 (u(t))Vu(i) G L 2 (Q). 

Since, due to (|2.15p . we know that Vu(t) G L 6 (fl) for almost all t, it is sufficient to check that 
6 _1 (u(f)) G L 3 (f2). Actually, this follows immediately from the proved fact that /(it) G L 2 (Q) (see 
the proof of the energy equality) and condition (|5.ip . 

Thus, we have verified that, for almost all t, m(u(t)) G iJ 1 (fl) and b(u)Vw(t) G L 2 (il). This 
justifies the first two equalities in f|5 -28[) . Note that the last one is obvious since w(t) G L 2 {yi) and 
Au(t) G L 2 (VL) for almost all t. Thus, we have verified that 

A(M(«(t)) J l) = Kt),A«(t)) 

for almost all t. Inserting the expression for w(t) from (11.21) into this identity, we end up with the 
desired entropy equality (|5.4p . Theorem 15-11 is proved. I 

Our next task here is to obtain stronger results on the attraction to the above constructed trajectory 
attr actor under the additional assumptions of Theorem l5.ll We start by stating a couple of corollaries 
that improve the results of Theorems 12 . 5 1 and 14.11 and simplify the construction of the trajectory phase 
space based on the energy and entropy inequalities obtained above. 

Corollary 5.4. Let the assumptions of Theorem 15.11 hold. Then, every weak solution of problem 
(jl.ip - (ll.2p satishes the dissipative estimate (|4.14p and, therefore, every weak solution automatically 
satishcs (|4. 16[) with C u = Q(£(u(0)). Thus, the condition (|4.16p in the definition of the trajectory 
phase space /C+ can be omitted and we may naturally consider /C+ just as the set of all weak solutions 
of problem (|l.ip - (11.2p . In addition, for every weak solution u, we have u G C([0, T], 77 1 (0)) and 

tt 1 -"ec([o,ri ) L 1 (n)). 

Indeed, we only need the entropy and energy (in) equalities in order to derive the dissipative estimate 
(|4.14p . Since these inequalities now hold for every weak solution, we have this estimate for every weak 
solution as well. The continuity properties stated in the corollary follow immediately from the energy 
equality. 
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Corollary 5.5. Let the assumptions of Theorem 1 5 . 1 1 hold. Then, the set B is bounded in JC+, in the 
sense of Definition 14.61 if and only if the set of the initial data {w(0), u £ B} is bounded in the energy- 
space Xfj_. 

We are now able to state our main result on the strong convergence to the trajectory attractor. 

Theorem 5.6. Let the assumptions of Theorem 15.11 hold. Then, the trajectory attractor A tr of 
problem (11.1[) - (|1.2[) is compact in C/ oc (R + , X) and the attraction property holds in that strong topology 
as well (remind that the compactness in C; oc (R+, X) means that the u-component of A tr is compact 
in Ci oc (R + , H 1 ^)) and the u 1-K -component is compact in Ci oc (R + , L 1 ^))). 

Proof. Let B C K + be the absorbing set introduced in the proof of Theorem 14.91 We claim that 
the set 

Bi := T(1)B 

is an absorbing set which is compact in the above mentioned topology (this is clearly enough for the 
proof of the theorem). Indeed, let {u n } C B be an arbitrary sequence of solutions. Then, since B 
is compact in 0™ eafe ; W e may assume without loss of generality that u n — > u in Q^ eafe ; where u also 
solves the problem (jl.ip ~ p.2p . To verify the above mentioned compactness, we need to check that 

u n -^u in C{[l,N],H l {ST)), u l - K -> u l ~ K in C([l, N], L 1 ^)), (5.29) 

for every N > 1. Furthermore, without loss of generality, we may check these convergences for N = 2 
only. 

To this end, we will use the proved energy equality which we will rewrite in the following form: 
i||Vu„(T)|| 2 + (F(u n (T)), 1) + (T(u n (T)), 1) - (g,u n (T)) 

T 

St\\d t u n (t)\\ 2 +t(b(u n (t))\7w n (t),Vw n (t)) dt 



T 1 

-||Vu„(t)|| 2 + (F(u n (t)),l) + (T(u n (t)),l) - (g,u n (t)) dt, (5.30) 

where T £ [1,2]. Our next task is to pass to the limit n /* oo in this inequality. First of all, thanks 
to the energy and entropy estimates and to the Aubin-Lions compactness theorem, we have 

u n ^u in C W ([1,2];V), (strongly) in C([l, 2]; H) n L 2 (l, 2; V), (5.31) 

and pointwise (a.e.). Thus, using the first convergence above and Fatou's Lemma, we see that 

h\\7u(T)\\ 2 < liminf \ \\ Vu n (T)f, (F(u(T)),l) < liminf (F(u n (T)), 1). (5.32) 

Next, thanks also to (|2T4]) - ([23)l . 

(T(«), 1) - (g, u) = lim (T(u n ), 1) - (g, u n ) strongly in C°([l, 2]). (5.33) 

Moreover, thanks to lower semicontinuity of norms w.r.t. weak convergence and to Ioffe's theorem, we 
also have 

St\\d t u(t)\\ 2 dt < liminf / St\\d t u n (t) || 2 dt, (5.34) 

n/-oo Jo 

t(b(u(t))Vw(t),Vw(t)) dt< liminf / t(b(uJt)Ww n (t),Vw n (t)) dt. (5.35) 

rv*oc J n 



Finally, in order to pass to the limit in (15 -30|) . we only need to prove that F(u n ) — > F(u) strongly in 
L 1 ([0, 2] x O). Actually, this follows from the uniform L 2 -bound of f(u n ), the pointwise convergence 
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u n — > u and the generalized Lebesgue theorem. Thus, we can take the supremum limit n oo in 
(|5.30p and obtain the inequality: 



lim sup T 



i||VMT)|| 2 + (FK(T)),l) 



< 



-T[(T(u(T)),l) - (g,u(T))] - / St\\d t u(t)\\ 2 - t(b(u(t))Vw(t),\7w(t)) dt 



-\\Vu(t)\\ 2 + (F(u(t)), 1) + (!>(*)), 1) - (g,u(t)) dt. 



(5.36) 



On the other hand, applying the energy equality in the form (|5.30l) directly to the limit solution u, 
and comparing with (|5.36[) . we obtain, for all T £ (1,2), 



limsupT -||V M „(T)|| 2 + (F(u n (T)), 1) 



< T 



-\\Vu(T)\\ 2 + (F(u(T)),l) 



whence, recalling (|5.32j) . we infer 

|K(T)|| V -> \\u(T)\\ v , \\F(u n (T))\\ LHn) -> \\F(u(T))\\ LHn) . 

This, together with the weak convergence u n (T) — > u(T) in V (cf. f|5.31[) ) and u n 
everywhere, implies the strong convergence 

u n (T)^u(T) in V, u 1 n ^ K {T)-^u 1 ~ K {T) in L^O), 



almost 



(5.37) 



for all T e [1,2]. This gives the strong convergence u n — > u in X pointwise in time. The desired 
uniform convergence (|5 . 29[) can be easily obtained using the standard contradiction arguments and 
applying the energy equality for u n (T n ) instead of u n (T). Theorem l5.6l is proved. I 

Corollary 5.7. Arguing in a similar way (and using also the entropy equality), one can verify the 
compactness and strong convergence to the trajectory attractor in all spaces involved in (|2.15[) - (|2.20j) . 



6 Separation from singularities and uniqueness 

In this section we prove that, in the viscous case 5 > 0, if k is large enough, then any weak solution 
becomes uniformly strictly positive for any t > 0. This is the object of the following 

Theorem 6.1. Let assumptions (|2.ip - (|2.6l) hold and let, additionally, S > and f3 — 0. In addition, 
let 

n>2s + 3ifd = 3, and k > s + 1 > 2 if d = 2. (6.1) 

Then, there exists a function Q; [0, oo) 2 —¥ [0, oo), monotone in each of its arguments, such that any 
weak solution u satisfies, for any e > 0, the separation property 

\\u- l (t)\\ L oo m <Q(E Q ,e- 1 ) fora.c.t>e. (6.2) 

The proof of the theorem will be given later in this section. As a consequence, we also have further 
time-regularization properties that imply uniqueness for strictly positive times as well: 

Theorem 6.2. Let the assumptions of Theorem 16.11 hold (in particular, let 5 > 0). Then, for any 
e > and any weak solution u there holds: 

w £ L 2 (e,T;V), (6.3) 
u e H 1 (e,T;V)nL°°(e,oo;H 2 (n)). (6.4) 

Moreover, in the class of weak solutions uniqueness holds at least for strictly positive times. 
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Remark 6.3. We point out that (|6.3)) - (|6.4D . which suffice to prove uniqueness, arc however not 
presumed to be optimal properties. Actually, thanks to (|6.2[) . (11 . If) is nondegenerate and (jl.2p is 
nonsingular for strictly positive times. Thus, by means of classical methods, one could easily prove 
that the solution u becomes, instantaneously in time, arbitrarily regular, provided of course that also 
the data 7 and g are smooth. 

As a consequence of uniqueness, we finally have 

Corollary 6.4. Let the assumptions of Theorem 16.11 hold (in particular, let 6 > 0). Then, the 
dynamical process generated by weak solutions admits the (strong) global attractor A in the standard 
sense (to be more precise, in the sense of semigroups with unique continuation). Namely, A is a 
compact and fully invariant subset of X^ such that, for any bounded set B C X^ there holds 

\im A x (u(t),A) =Q, (6.5) 

t yoo 

uniformly with respect to weak solutions u such that u(0) G B. 



6.1 Proof of Theorem 16.11 in the 3D-case 

We consider equation (|1.2p rewritten in the form (|5.19l) , where we assume 5 — 1 for simplicity. We also 
assume s > 0, the case s — being simpler since one can directly take advantage of the L 2 (0, T; L 6 (il))- 
regularity of w. Then, the proof is based on a suitable version of the Moser iteration argument, i.e., 
we will take v > 1 and test (|5.19l) by vz v ~ x for increasing exponents v. We have to remark that this 
procedure, apparently having a formal character since the above test function could grow very fast 
and, hence, have insufficient regularity, can be easily justified simply by truncating z at some level 
K and then letting K /*■ 00. In particular, the argument does not require any approximation of the 
equation and, hence, works for all weak solutions in the class introduced in Definition 12.31 

That said, testing (|5.19l) by vz v ~ x and integrating over Q v :— ft x (t u ,T), where the "initial" 
time t„ will be chosen later on, we then have 

J" + JJ Q z K+v+1 < \WuW L « {n) + "JJ q <^ +1 - v jJ wz " +1 ' ( 6 - 6 ) 

where we have set 



Adding now 



■K ■= \\4Ur„,T;L» m + \\VZ»I 2 \\1 W ,H)- (6-7) 

to both hands sides of (|6.6[) . in order to recover the full V^-norm of z v l 2 , and setting 

K '■= \\ Z WL^{T V ,T-L"(n)) + c n\\ z \\Lv{T, y ,T;L 3 »(n)) ^ IMIiWstQ,,): (6-9) 

where cq is a suitable embedding constant, we then arrive at 

z K+ » +1 < \\z(r v )\\l„ m + uC\\z\\l+l HQv) + \\z\\l„ {Qv) v JJ^ wz» +1 , (6.10) 

and we have to provide a bound for the last term on the right hand side. 

To do this, let (p,p*) and (q, q*) be two couples of conjugate exponents with p < 2. Then, 
using also (|5.22l) . we obtain 



wz v+1 < v\\w\\ LP[T ^ T . Lq[n)) \\z v+1 \\ LP * {TvT . Lq * m 

< v{\\w\\LP{T»,T;V{n)) + \\Vw\\lp{Q v ))\\z V+1 \\ L v* ( Tv ,T;Li* (n)) 

< v(C+ IIVwIU,^))^^ 1 !!^.^^,.^, (6.11) 
provided that we choose q — 3p/(3 — p), so that TV 1,P (S1) C L q (Q) continuously. 
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Then, the term with Vw is estimated this way: 

l|VHU*(Q,) = \\u s/2 z s/2 Vw\\ LP{Qv) < Wu^VwU^qJIz^W^^ 

<C||z s/2 || i , (6.12) 

thanks also to (|5.2ip . Thus, collecting (|6.11l) and (|6.12j) . we have 

- ' 1 < Cv(l + ll*"\^ w J \\4^^ M ^ my (6-13) 



wz 



Now, let us assume to know a bound of the term I n -i = I Vn _ x from the preceding step of the iteration. 
Then, thanks to the last inequality in (|6.9[) , we can use it to estimate the term in brackets so to have 

~ V ffq ^^( 1 + ^-0ll a ll^ 1 ) (rv , r;i «-c« + i)( n)) . ( 6 - 14 ) 

provided that one chooses p as follows: 

p 5i/ n _i 1 1 3s 

i.e., - = - + — . (6.15) 



2 — p 3s pi 10i/ n 

This gives in turn 

1 1 3s 1 1 3s 1 5 3s 



q 6 10v n -i P* 2 10z/„-i 9* 6 10z/ n -i 



(6.16) 



Then, we have to take v — v n in a way suitable for the next step of the iteration. The choice is 
dictated by the exponents of the last term in (|6.14j) ; namely, v = v n should be close enough to v n —\ 
in order that term be still controlled by l n -\. Using interpolation, we require that, for some 6 G [0, 1], 

11-'-' 1 ^ ' (6 . 17) 



p*(v n + l) oo v n -\ q*(u n + l) v n -i 3^„_i 
To compute 9, we first take the quotient of the above equalities and then use f|6. 16|) . This gives 

3-26* p* 25j/„_i - 9s 



36* q* 15f„_!-9s' 

whence 



(6.18) 



15f„_i-9s v n -i , 7i/ n _i - 3s - 6 / c1ri \ 
g= 35^-15, ' ^ ^ = W" 1 = 6 ' ( ° 9) 

where the second of (|6.16l) has also been used. 

Thus, it turns out that v n > v n -i provided that v n ~i > 3(s + 2). Thus, in order the above 
iteration could be performed, we need to find some n$N and some Vn > 3(s + 2) such that, for any 
e € (0, 1), there holds 



{\\4T^T;L^ m + ^W4T^ { e,T;L^ m ) " < ( 6 - 20 ) 



where Q is a computable monotone function (whose expression can depend on the magnitude of the 
initial data and of T). 

Let us pospone the verification of (16.201) and let us now see that, for n > n, the induction 
principle can be applied. Coming back to (|6.10p . we then have 

K n < IKOII£n(n) + Cv n {l + \\z\\lltl {Qn) ) + Cv n l£? + \ (6.21) 

where we wrote n in place of v n in some subscripts and assumed w.l.o.g. I n -i > 1- Thus, extracting 
the f n -th root and noting that v n + 1 < 5^„_i/3, we obtain 

In < Mr n )\\ L , Hn) + (C Vn )^rt_ 1: where Vn := 2l/ " + s + 2 ( 6 . 22 ) 
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and C is independent of n. Moreover, for (arbitrarily small) e G (0,1), given r„_i we can choose 
T n G [t„-i, t„_i + en^ 2 ] such that 

5"r»-l B"n-1 „2 /.T„_i+£ t l- 2 5^„_! 2 5v„_! 

lk(r B )|| L ^ (n) < 4<r n )\\ <c— / Ml .5^ <c-/ n 4 . (6.23) 

v ' l — 3 — (n) e Jr„_i l — s — (n) 6 

Thus, (|6.22l) can be rewritten as 



,2 



c- r*- 1 +(CO- (6-24) 



4 < 

whence a standard computation permits to pass to the limit w.r.t. n /* oo. Since lim n ,/* 0o r„ exists 
and is less or equal than ce, we then obtain 

IMU~(Qx(e,T)) <Q(e- x ), (6.25) 

for Q as in (|6.20[) . as desired. 

Thus, to conclude the proof it only remains to check that (|6.20[) holds. To do this, we come 
back to (|6.10l) and use now the i K+!/+1 -norm to estimate the right hand side. Proceeding as above, 
we still arrive at (|6.13p , where now we have to take 

P K + fn-l + 1.11/ s \ 

i.e., - = - l + — — . (6.26) 



2— p S p 2\ K + Vn-i + 1 

Thus, we obtain 

1 1/1 s \ 11/ s ^ 11/5 s \ 

q ~ 2 V3 + k + v n -i + 1/' p* ~ 2 V ~ K + v n -i + 1/' ~~ 2 V3 ~ k + z/„_i + 1/' ' ' ' 

whence we get the analogue of (|6.14[) . i.e., 

<CKl + Al\)||z||^ ( „ +1)(T ^ T;ig , ( „ +1)(0)) , (6.28) 
where 

:= IMI£~ (T „,x ; ^(n)) + l|V^ /2 ||! 2(T „, T;ff) + IMl£?+ + V) (6-29) 

and A„ := A„ n , as before. Then, we still have to choose v = v n in a suitable way. Similarly as before, 
we require that for some 8 G [0,1] it is 

1 i a a 1 i a q 

(6.30) 



P*(l>n + 1) OO K + V n -i + 1 g*(z/„ + l) I/ n _! K + I/„_i + l 

To compute 6>, we first take the quotient of the above equalities and then use (|6.27[) . This gives 

0v n -i + (1 - 6){Vn-l + K + 1) _ p* _ 5(l/ n _i + K + 1) - 3s 



9v n -i q* 3(z/„_i + k + 1 — s) 

whence 



(6.31) 



and, from the first of f|6 . 30[) . 



0~ 1+ 3(i/ n _i + K + l)-3a (6 ' 32) 



u„ - 1 + k + 1 5 1 , , , 

1 = + ~(k-s-1), (6.33) 



Op* 6 2 



whence it is clear that v n > is n -i if and only if v n -\ < 3(k — s — 1). Then, proceeding similarly 
with the previous part of the iteration, if we start knowing a bound of h. Va for some i/q > 1, then we 
can reach, in a finite number n of steps, any v-^ < 3(k — s — 1). Since we also need > 3(s + 2) 
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from before, this leads to the compatibility condition 3(s + 2) < 3(k — s — 1), that is equivalent to 
assumption (|6.1[) . 

Thus, the proof is concluded provided that we find vq > 1 to start the argument. Actually, 
we can test (|5.19p by z L for small i > 0. We obtain, for Q = SI x (0, T), 

Jltl +JJ q - K+2+t < INo||itt ( a) +JJ ^ 2+ ' -JJ ^ 2+L (6-34) 
and, being n > 3 by (|6.1|) . it is clear that, at least for i < 1, 

>z 2+t - // ^z 2 +' < 1 // z K+2 +' + C, (6.35) 

thanks also to (|5.2[) . Hence, we can take z^o = 1 + i for arbitrary i £ (0, 1), which concludes the proof 
in the case d = 3. 

6.2 Proof of Theorem 16.11 in the 2D-case 

The proof is carried out by the very same scheme used in the 3D case, the differences being limited to 
the exponents related to use of interpolation and embeddings. Thus, we limit ourselves to point out 
these differences. Now, in place of (|6.9p . we have 

K ■= \\4U M L^n)) + \\^\\h(r v ,T;V) > \Wv»{Q„y (6-36) 

Thus, taking p £ (1, 2), we have q = 2p/(2 — p), so that, to control the right hand side of (|6.14p . we 
need to choose p so that 



p 2v n - X .11- 

Le -> Z = o [ 1 + ^7. > ( 6 - 37 ) 



2 — p s p 2 V 2v. 

whence we obtain 

1 s 11/ s \ I 



\ = l(i-r L -\ ^ = 1 -j ± - (6-38) 

n* 2 V 2v n -^ / q* 4^„_i 



q 4l^„_i' p* 2\ 2v n - 1 ) 1 q* ->r„ 

and, correspondingly, 

1 1-9 9 1 1 



Thus, 



whence 



p*{v n + l) oo 2^„_i' q*{v n + l) v n -i 2v n ^i 

2 — 9 p* 4z/„_i — s 
9 q* 2v n -i-s' 



(6.39) 
(6.40) 



2v n _i-s 3 s + 2 

9 = - , and v n = -v n -\ — , (6.41) 

3*Ai-i — s 2 2 

so that we need to find Ufi > s + 2 in order the procedure works. 

To do this, we proceed again as before and, choosing p as in (|6.26l) . the other exponents are 

then given by 

1 s 11/ s \ 1 s 

q ~ 2(k + v n -i + 1) ' p* ~ 2\ ~ k + v n -i + 1/' <F ~~ ~ 2(k + i/„_i + 1)' ^ ' ^ 
Then, taking e [0, 1] as in (I6.30[) . we now arrive at 

9u n _ 1 + (I - 9)(v n _ 1 + k + I) _P^_ 2(i/„_i + « + !) -s 

6v n -i q* Vn-x+K+l-S 

whence 

- = 1H , and v n = v n -i + -(k - 1 - s), (6.44) 

" ^71-1 + K + 1 — S I 

so that it is v n > v n -\ if and only if k > s + 1, i.e., (|6.1[) holds. Thus, we can arrive in some finite 
number n of steps to have > s + 2 provided that we can start as before from vq = 1 + 1 for some 
(small) l > 0. Actually, we can now take t = k — 2, which is strictly positive thanks to (16.11) . Thus, 
(|6.34l) can be repeated without any variation and, of course, we still have (|6.35[) thanks to Holder's 
and Young's inequalities. The proof is complete. 



.43) 
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6.3 Proof of Theorem 16.21 and Corollary 16.41 

Again, we just consider the case d = 3, the case d = 2 being simpler. First of all, we deduce further 
regularity of weak solutions. Actually, thanks to (16.21) . u is uniformly separated from for any time 
t > e > 0, e>0 being arbitrary. Then, (jl.ip becomes in fact nondegenerate and the energy estimate 
gives the improved regularity (|6.3I) . Moreover, the term f(u) in (11.21) is now smooth and we can 
apply the linear parabolic theory (or test (|1.2p by —(t — e)Ait t and perform standard computations) 
to deduce (jB~4|) . 

At this point, rewriting (jl.ip as a family of time-dependent elliptic problems, namely 

- Aw = — !— ( - u t + b'(u)Vu ■ Vw) , (6.45) 
b{u) v ' 

relations (|6.4[) and (|6.3p permit to see that the right hand side belongs to L 2 (e, T; L 3 / 2 (fJ)), whence 
we obtain 

w € L 2 (e, T; W 2 > 3/2 (Q)) C L 2 {e, T; W 1 ' 3 ^)). (6.46) 

To prove uniqueness, we can now consider a couple of solutions u%, u 2 , set u :— u\ — u 2 (and, 
correspondingly, w :— w\— w 2 ) and take the difference of equations (ll.ll) - (jl.2p to obtain 

u t - div(6(ui)Viw) = div ((6(ui) - b(u 2 )) Vw 2 ) , (6.47) 
w = Su t -Au + W'{ui) - W'(u 2 ), (6.48) 

where W = f + 7 can be thought to be globally Lipschitz in view of the strict positivity of u\ and 
U2- Then, we test (|6.47[) by w and (|6.48p by ut- We obtain, for some c, a > 0, 

^HV M || 2 + a||VH| 2 + <% t || 2 

< c||u||||«t||+ / |(6(«i)-6(«2))V«;-Vtfi2| (6-49) 
Jn 

and we can estimate the last term as follows: 

|(6(«i) - b(u 2 ))Vw ■ Vw 2 \ < ^\\Vw\\ 2 + c\\u\\ 2 L6{n) \\S7w 2 \\ 2 L3{n) 

< fl|V W || 2 + C || U || 2 y ||V W2 || 2 3(a) . (6.50) 

Thanks to (|6.46p . we can then apply Gronwall's Lemma to (I6.49p . which gives the assert. At this 
point, Corollary E3] is an immediate consequence of the uniqueness property and of the general theory 
of infinite-dimensional dynamical systems [2 122j . 
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